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I PACS. 71.10.Fd - Lattice fermion models (Hubbard model, etc.). 

■ PACS. 71.30.+h - Metal-insulator transitions and other electronic transitions. 
PACS. 71.10.Hf - Non- Fermi-liquid ground states, electron phase diagrams and phase transi- 
tions in model systems. 
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Abstract. - The Tomonaga-Luttinger parameter Kp determines the critical behavior in quasi 
one-dimensional correlated electron systems, e.g., the exponent a for the density of states 
near the Fermi energy. We use the numerical density-matrix renormalization group method 
to calculate Kp from the slope of the density-density correlation function in momentum space 
at zero wave vector. We check the accuracy of our new approach against exact results for 
the Hubbard and XXZ Heisenberg models. We determine Kp in the phase diagram of the 
extended Hubbard model at quarter filling, tIc = 1/2, and confirm the bosonization results 
Kp = = 1/4 on the critical line and K'^^^ — n^/2 — 1/8 at infinitesimal doping of the 
charge-density-wave (CDW) insulator for all interaction strengths. The doped CDW insulator 
exhibits exponents a > 1 only for small doping and strong correlations. 
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' Quasi one-dimensional correlated electronic systems offer the unique opportunity of a de- 

tailed comparison between theory and experiment. Some models for correlated lattice electrons 
can be solved exactly, e.g., the XXZ Heisenberg chain and the Hubbard model, and many of 
the physical properties of Hubbard-type models are known from analytical and numerical 

■ studies. Moreover, their low-energy properties turn out to be universal, i.e., they belong to 
^ . the generic class of Tomonaga-Luttinger liquids (TLL). It has been a challenge to detect the 

I ' signatures of TLL in experiments on metallic single- walled carbon nanotubes [1] , anisotropic 
^ ! transition-metal oxides like PrBa2Cu408 [2], Bechgaard salts Uke (TMTSF)2C104 [3], and 
g , TTF-TCNQ [4]. 

^ , The low-energy properties of TLL are characterized by few quantities, most notably the 

^ ' TL parameter Kp which determines the long-range behavior of the density-density correlation 
, function and the exponent a for the density of states D{lu) near the Fermi energy [5], 

D{uj^O)^\ujr , a={Kp + K;^-2)/4. (1) 

A value a > 1 was estimated for TTF-TCNQ [4] and in early works on (TMTSF)C104 [3] 
which would imply Kp < 3 — 2^/2 « 0.17. More recent measurements report Kp = 0.23 for 
(TMTSF)C104 [6] and Kp = 0.24 for PrBa2Cu408 [2]. In the Hubbard model, Kf > 1/2 so 
that < 1/8. Therefore, if electronic correlations are indeed responsible for the algebraic 
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behavior of the density of states near the Fermi energy, a > 1/8 is only possible when the 
long-range parts of the Coulomb interaction are taken seriously. Unfortunately, the calculation 
of Kp for corrclatcd-electron models is very difficult. 

In this work we employ the density-matrix renormalization group (DMRG) method for a 
reliable numerical calculation of the TL parameter. In order to demonstrate our approach 
we consider interacting spin-1/2 electrons on a chain with an even number L sites. The 
electron density is n| = — n/2 = N/(2L). In the absence of a Peierls dimerization, the 
extended Hubbard model provides an appropriate description for correlated electron systems 
in one dimension, 

H =-t^ (c;+^c;+i,<, + c^i,<,c;,a) ^;,T"i,i + V^{ni-n) (n;+i - n) , (2) 

l,a I I 

where c/"^ {ci,a) is the creation (annihilation) operator of an electron with spin a on site /, 

'^;,o- = is the number operator, and = n; | + n; |. Moreover, t is the electron 

transfer integral between neighboring sites, U is the strength of the Hubbard interaction, 
and the nearest-neighbor interaction V models the long-range part of the electrons' Coulomb 
repulsion. 

The density-density correlation function is defined by the ground-state expectation value 

1 ^ 

CNN(^) = - J^if^l+rfll) - {m+r){ni} , (3) 
1=1 

where hi = J2a ^i,^- counts the electrons on site I. We have C^^{r) = C^^{—r) due to 

inversion symmetry, and periodic boundary conditions apply. 

Using conformal field theory it can be shown [7, 8] that the asymptotic behavior for 1 <C 
r <C i is given by 

where fcp = n7r/2 is the Fermi wave number, and A is a constant. For spinless Fermions, the 
first term should be multiplied by 1/2. Field theory further predicts [9] Kp{Vc{U); ric) = n1 for 
an interaction-driven metal-insulator transition. For example, we have Kp{yc{U); 1/2) = 1/4 
for the CDW transition in the quarter-filled extended Hubbard model. In contrast, for the 
density-driven metal-insulator transition the value from field theory is KpiV > Vc(J7),n~) = 
nl/2, e.g., for the extended Hubbard model field theory predicts Kp{V > Vc{U), 1/2") = 1/8 
in the infinitesimally doped CDW phase at quarter filling. 

As seen from the exactly solvable cases discussed below, we expect that Kp{U, V; ric — S) 
can be expanded in a Taylor series in the doping 6 = ric — n <^ 1, 

i.,(f/,y,n.-.) = | + ^ + .... (5) 

Moreover, since the critical line is a line of Kosterlitz-Thouless transitions, we expect that 
the prefactor diverges exponentially in the vicinity of the critical line, 

My^Km~exp(^^,-|^) . (6) 



Satoshi Ejima et al.: Tomonaga-Luttinger parameters 



3 



Therefore, close to the critical line the value Kp — of the CDW insulator cannot be 

observed for finite doping. 

Field theory does not quantify the convergence radius of (O or the region in which (0 holds. 
To answer this question, the TL parameter must be calculated for the extended Hubbard model 
with accurate numerical methods. We can extract Kp via Fourier transformation, 

L 

CNN(g) =^e-""-CN^(r) , (7) 

r=l 

with < g < 2ti. By construction, C^^{q = 0) = 0. For the derivative at g = one finds in 
the thermodynamic limit [10,11] 

Kp^^ lim ^l^^lM . (8) 

9-.0+ q ^ ' 

In numerical simulations we treat finite systems. There, Eq. ijSJ translates into 

Kp = hm f ^\ , 



Several groups have calculated the density-density correlation function in position space for 
Hubbard-type models. After Fourier transformation they obtained Kp from Q, see, e.g., 
[10-13]. The main problem of this approach lies in the accurate calculation of C^^(r) from Q 
for large distances. The accuracy of the correlation function becomes worse as the distance r 
increases which severely limits the precision of the Fourier transform C^^{q), especially for 
small q. ^ 

In this work we calculate C^^{2tt/L) directly in momentum space. To this end we define 

N[q) = j{^o\n{<l)ni-<l)\'fo) (10) 

for q = 27rm/i (m > 1), where n(q) is given by 

niq) = n+i-q) = ^ ^""'^'"'^^^^+,'5,. . (11) 
La 



Here, Vc ^ {L + l)/2 denotes the central position of the chain. Note that N{q) and C^^{q) 
are different. It is only in the thermodynamic limit, when boundary effects are absent, that 
they become identical. Therefore, 

L _ /2tt 



Kp{L) = -Ni^—j , Kp^ Ihn^KpiL) . (12) 

The important idea is to target not only the ground state in the DMRG procedure but also 
the state j^'g) ~ n{—q)\'i/Q). In this way, a precise DMRG calculation of N{q) and of Kp 
from (|12|l becomes possible. 

We illustrate the accuracy of our method for the Hubbard model [V = in |5J], for which 
Kp{U ; n) is known from the Bethe Ansatz solution [8]. We investigate systems with L < 128 
and open boundary conditions. The number of density matrix states kept is m = 1500, so 
that the maximum truncation error is 3 x 10~^. 
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Figure 1 - Left: Kp{L) as a function of the inverse system size in the one-dimensional Hubbard 
model ai n — 0.1 for U/t—1 (circles), U/t — 2 (squares), U/t — 6 (diamonds), U/t = 10 (upward 
triangles), and U/t — oo (downward triangles). Diamonds are exact values from the Bethe Ansatz in 
the thermodynamic limit. Lines are 4th-order polynomial fits. 

Right: TL parameter Kp in the one-dimensional Hubbard model as a function of the density for for 
U/t — 1,2,6, 10 (from top to bottom). The full lines are exact values from the Bethe Ansatz, stars 
mark the results from DMRG. 

In Fig. we show Kp{L) as a function of the inverse system size for several values 
of U/t. The band filling is fixed at n = 0.1, which, apart from the hmit n ^ 1, is the most 
difiicult parameter region in this model because Kp changes significantly as a function of the 
interaction strength. For all values C/ > 0, Kp{L) is found to decrease monotonically as a 
function of inverse system size, so that we can extrapolate Kp to the thermodynamic limit 
systematically by performing a least squares fit of Kp to a polynomial in 1/L. 

In Fig. ^ we compare our results for if^^^*^ with those from Bethe Ansatz for various 
filhngs and interaction strengths. The relative error \Kf^^'^ - iv:cxact|^^exact jg ^^^^^ q 
for all DMRG data shown. We reproduce the exact results with a much better accuracy 
than exact diagonalization [8], the DMRG method [11] based on the Fourier transformation 
formula ©, and the calculation of Kp from the compressibility and the charge velocity [14]. 

For the Mott-Hubbard transition at Uc = 0+ the Taylor expansion © applies [15] {uc — 1, 
f{U) = h{U.t))) 

Fig. m shows our results for the infinitesimally doped Mott-Hubbard insulator, rih = 2/L, 
for various values of the interaction strength. The approach is valid because we calculate 
N[2-k/L) so that the momentum q = 27r/L transferred from the ground state |\l/o) to \^q) is 
of the same order of magnitude as the phase shifts induced by the introduction of two holes. 
For U ^ 2t our systems are larger than the correlation length of the system, and we reproduce 
Kp{U > 0, n ^ 1) = 1/2 numerically. Finite-size effects become prominent for U ^ 0, and it 
is difficult to recover the exact value from numerical calculations in that limit. 

As our second test we study the extended Hubbard model at quarter filling, n = 1/2, 
for U ~ 00, which can be mapped onto the exactly solvable Heisenberg XXZ chain. The 
parameter Kp for the Hubbard model at L/ = cx) from the Bethe Ansatz is {V <Vc = 2t) 



4arccos[-F/(2t)] 



(14) 
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Figure 2 - Kp{L) as a function of the inverse system size for the infinitesimally doped Mott-Hubbard 
insulator, n — 1 — 2/L, for various interaction strengths. It extrapolates to the exact value Kp = 1/2 
for all U/t. 



In Fig. we show Kp as a function of y for C/ = cx) together with the exact resuh. For this 
system we use periodic boundary conditions because m = 2000 density-matrix eigenstates 
are enough to calculate Kp[L) with high precision, and finite-size effects are much smaller 
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Figure 3 - Left: TL parameter Kp for the extended Hubbard model at U/t — oo, V < 2t from 
DMRG (stars) and exact Bethe Ansatz (full line). Inset: Kp{L) as a function of inverse system size 
for periodic boundary conditions. 

Right: Contour map for the TL parameter Kp in the U-V plane of the extended Hubbard model 
at quarter filling. The bold line represents the boundary of the metal-insulator transition. The 
infinitesimally doped charge-density- wave insulator (CDW) has Kp — 1/8. The shaded area indicates 
the region with an exponentially small gap. 
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for periodic than for open boundary conditions, see the inset in Fig. Relative errors 
I^DMRG _ ^oxact|/^cxact g^^^ ^^i^^ o.5%, evcn for V = 1.95t where the system is close to the 
CDW insulator. Again, the predictions (O, ((SJ of field theory apply. Here, the Bethe Ansatz 
solution [16,17] gives (ric = 1/2, 2g{V) = h{oo,V), cosh(7) = V/Vc) 

fy. ^ l + 2E,TLi(-l)"/cosh(n7) 
' l + 2E,T=i[l-tanh(n7)] ^ 

giV^V^) = -^expf , ] 

As an application, we study the extended Hubbard model at quarter filling. In Fig. (Sja, 
we show the phase diagram together with the contour lines for the TL parameter Kp. Three 
different phases are found, namely, a 'superconducting phase' {Kp > 1), a metallic phase 
(1/4 < Kp < 1), and a 4fci?-CDW insulator beyond the critical line. The results are in good 
agreement with previous works [18, 19]. 

On the CDW transition line we find Kp = 1/4, and Kp^"^ = 1/8 for the infinitesimally 
doped CDW insulator, in agreement with field theory [9]. The parameter region where finite- 
size effects are prominent due to an exponentially small gap is shaded in Fig. Outside this 
region we can determine the TL parameter reliably. 



(15) 
(16) 




Figure 4 - TL parameter Kp as a function of V/t for U/t = 6, 10, oo at filling n — 11/24 (from top 
to bottom). The full line is the result H15^ . the dashed horizontal line marks Kp = 3 — 2y/2 ~ 0.17. 
Inset: U/t = oo for fillings n = 21/48, 22/48, 23/48 (from top to bottom). 

The results for S — 1/24 w 4% are shown in Fig. 0] as a function of V/t for U/t = 6, 10, oo. 
Deep in the CDW phase {U/t > 5, V/t ^ 6) neither U nor V have a large influence on Kp. 
Fig. 0] also shows, however, that Kp < 0.17 and thus a > 1 requires very large interaction 
strengths, or very small doping, 5 ^ 2%. 

In conclusion, we developed an accurate numerical DMRG method to obtain the TL pa- 
rameter Kp for Hubbard-type models. We demonstrated its accuracy for the Hubbard model 
and the XXZ Heisenberg chain (spinless fermions). We presented an accurate phase diagram 
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of the t-U-V model at quarter filling and verified the field-theoretical predictions for Kp. We 
also showed that a critical exponent a > 1 is only possible for a lightly doped CDW insulator 
with a sizable gap. 

* * * 
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